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0.  Introduction 


KIDDER  (1976)  studied  the  homogeneous  (i.e.,  uniform) 
isentropic  compression  of  hollow  shells  imploded  by  the 
conversion  of  absorbed  laser  energy  into  mechanical  work. 

(The  characteristic  feature  of  uniform  motion  is  that  the 
radial  velocity  at  a  given  Instant  of  time  is  everywhere 
proportional  to  the  radial  distance  from  the  center  of 
symmetry.)  Using  this  model  he  was  able  to  estimate  the 
optical  power  required  to  achieve  a  given  measure  of  Inertial 
confinement  for  a  hollow  laser-driven  pellet.  The  calcul¬ 
ations  showed  that  using  a  large-aspect-ratio  target  (radius 
much  greater  than  thickness)  reduced  the  power  requirements 
compared  with  those  for  a  solid  spherical  target  by  as 
much  as  a  factor  of  five.  The  results  were  strongly 
influenced  by  the  possibility  of  disruption  caused  by 
Rayleigh-Taylor  instabilities,  which  imposed  an  upper  limit 
on  the  aspect  ratios  which  could  be  employed. 

The  time-dependent  basic  state  in  the  model,  a  type  of 

self-similar  motion,  is  of  considerable  theoretical  interest. 

Although  simple,  it  subsumes  a  number  of  the  principal 

physical  features  of  the  problem,  and  represents  a  nontrivial 

solution  of  the  (nonlinear)  equations  of  motion.  Such 

solutions  appear  to  have  a  range  of  applicability  well 

beyond  what  one  might  have  expected,  apparently  because  the 

change  of  scales  implicit  in  strong  compression  is  large 

enough  so  that  the  dependence  on  variables  other  than  the 
Note:  Manuicript  submitted  October  17, 1979. 


1 


similarity  variable  becomes  weak.  Similarity  solutions 
of  ’^ais  and  other  types  have  accordingly  been  found  useful 
in  a  wide  range  of  studies  of  laser-driven  target  behavior 
[see,  e.g.,  A»ISIMOV  and  INOGAMOV  (1974)]. 

It  is  noteworthy  that  KIOOER  (1976)  found  it  possible 
,'.o  carry  out  in  closed  form  the  analysis  of  the  linear 
stability  of  Incompressible  irrotational  perturbations  about 
the  self-similar  basic  state.  Thus  he  was  able  to  determine 
the  evolution  of  the  Ray leigh-Tay lor  instability  exactly, 
without  recourse  to  numerical  approximation,  even  though  the 
perturbed  equations  of  motion  contain  coefficients  which 
depend  on  position  and  time,  and  the  perturbations  are  there¬ 
fore  not  the  familiar  exponential  functions  of  time. 

Somewhat  later,  the  present  authors  independently  dev¬ 
eloped  in  a  different  application  (BOOK  and  BERNSTEIN,  1978) 
a  technique  for  determining  the  evolution  of  all  the  modes 
(including  compressible  and  rotational  ones) .  Now  we  wish 
to  apply  this  technique  in  order  to  generalize  KIDDER'S 
(1976)  results  in  three  respects:  (i)  arbitary  choice  of 
Y,  the  ratio  of  specific  heats:  (ii)  nonuniform  entropy; 
and  (ill)  general  perturbations. 

As  discussed  below,  we  obtain  the  following  principal 
results:  (i)  The  Rayleigh-Tay lor  modes  become  more  unstable 
as  Y  decreases  (this  is  true  of  all  modes) ;  (ii)  they  are 
unaffected  by  changes  in  the  shell  density  profile; 
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(ili)  compressible  modes,  which  are  convective  in  nature, 
are  destabilized  by  the  occurrence  of  regions  of  "super- 
adiabaticity"  [d(pp”^)/dr  >  O]  ,  but  grow  less  rapidly 
than  incompressible  modes,  particularly  at  |.nfinitely  large 
mode  number,  where  the  Rayleigh-Taylor  growth  rates  diverge. 

The  paper  is  organized  as  follows:  in  SI,  we  write  down 
and  solve  the  equations  describing  the  unperturbed  motion. 

In  S2  we  derive  the  equations  describing  linear  perturbations 
about  the  basic  state.  The  time  and  space  dependence 
separate.  The  time-dependent  part  is  solved  analytically  in 
S3  for  general  values  of  the  separation  constant.  In  S4 
we  solve  the  spatial  part  of  the  problem  for  incompressible 
irrotational  modes.  In  S5,  utilizing  a  variational  principle, 
we  solve  the  spatial  equations  for  the  case  of  compressible, 
rotational  modes.  In  S6  we  illustrate  the  results  so  obtained 
with  a  number  of  examples  which  admit  of  analytic  solutions 
of  the  spatial  problem.  (The  variational  principle  of  S5 
is,  however,  well-suited  for  numerical  determination  of  the 
separation  constant) .  We  conclude  in  §7  with  a  discussion 
of  our  results  and  some  applications  to  targets  having  typical 
parameters . 


1.  The  model 

We  start  with  the  equations  of  ideal  hydrodynamics  in 
the  form 

p+v7p-0,  (1. la) 


I 

I 

! 
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(1.6) 


^Pq  -2 

-  •’o'  ^  • 

Here  a  constant  t  with  units  of  time  results  from  separating 
the  spatial  behavior  (1.6)  and  the  equation  determining 
f  (t)  : 


(1.7) 


The  sign  in  (1.7)  is  chosen  to  yield  negative  acceleration 
(implosion).  As  a  consequence  of  this  choice,  Pq(>^)  incre¬ 
ases  monotonically  outward.  From  (1.5)  it  follows  that  the 
pressure  applied  to  the  outer  boundary  of  the  shell  at 

r  ■  JfQ+d  must  grow  in  time  proportionally  to  f 

general  a  quadrature  can  be  carried  out  on  (1.7),  yielding 


2  *2 

T  f  ■  -  2071  f 


(1.8a) 


for  Y~l>  And 

'f2  -  ^  (f"®-l)  (1.8b) 

a 

otherwise,  where  a  «  3(y-1).  For  y  ■  5/3,  (1.8b)  can  be 

solved  to  give 

f (t)  -  (1-t^/T^)'  .  (1.9) 

For  other  values  of  y  the  solution  is  easily  found  numer¬ 
ically.  In  all  cases  f(t)  is  approximately  parabolic  near 
t  ■  0,  and  vanishes  with  infinite  slope  at  some  time  t^  of 
order  t . 

Although  the  form  of  the  time-dependent  function  f(t) 
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is  determined  by  specifying  y,  there  is  complete  freedom 


in  the  choice  of  nonnegative  PQ(t)  (jor  equivalently,  p^lr)  , 
provided  the  latter  satisfies  Pq(z')  >  0]*  Por  example,  if 


(1.10a) 


where  H  is  an  arbitary  constant,  then 


/  ,  -  /  2  ,  2  , ,  *M-1 

Po(r)  -PQ{r  Aq  -1) 


(1.10b) 


«  M  2  2 

with  Po"^0^0  /2(^+1)t  .  The  case  h  »  0  corresponds  to  a 
shell  of  uniform  density  and  quadratic  pressure.  When 

entropy  s(r)  becomes  constant.  This  is  the 
model  studied  by  KIDDER  (1976)  .  A  particularly  simple 
limiting  form  of  (1.10),  which  will  be  employed  in  §  6,  is 


/X  2  ,  2.  H 

PQ(r)  -  p(r  /r^  )  » 


(1.11a) 


,  ,  2,  2.M-1 

PQ(r)  -  p(r  /r^  ) 


(1.11b) 


Note  that  in  systems  described  by  this  model,  all 
motion  is  inward  —  there  is  no  ablation  region  or  blowoff, 
nor  does  heat  conduction  play  a  role.  Hence  application 
of  the  model  is  restricted  to  the  cold  inner  regions  of  a 
shell  or  pellet. 


2.  Linear  perturbations 


We  follow  BERMSTEIM  and  BOOK  (1976)  and  BOOK  (1978, 
1979}  In  deriving  linearized  equations  written  in  Lagranglan 
variables  for  the  development  of  a  small  perturbation  about 
the  solutions  of  the  basic  equations.  As  a  result  of  the 
perturbation,  an  element  of  fluid  whose  unperturbed  motion 
was  given  by  ^  (j»t)  is  located  atJR(j:,t)  +jJ(jc,t).  The 
first-order  displacement  satisfies  the  linearized  form 


of  (1.1b), 


Pi  +  Pl5-  -  ^rP, 


(2.1) 


where  the  subscript  1  denotes  first-order  quantities. 


Substituting  for  p  and  p.  from 

X  * 


P  *  -  C 

1  R  »» 


(2.2) 


Pi  -  tt^l  -  -  i. 


(2.3) 


and  using  (1.3)  to  replace  the  gradient  with  respect  to 
JR  by  the  gradient  with  respect  to  according  to  f  ^ 


we  find 


tVY-1  ^  j  7.^  -  t2yPq’^7(Pj,7.0  +  7J.r.  (2.4) 


In  (2.4)  we  have  omitted  the  subscript^  on  V.  For 
convenience  in  what  follows  we  will  also  omit  the  zero 
subscript  from  and  p^. 

Iiihe  the  unperturbed  equation,  (2.4)  is  separable. 


Writing 
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-  ^(r)  T(t) 


(2.5) 


wa  have 

^  +  T^YP"^  V  (pV.^)  +  {y^) -  0  (2.6) 

and 

,2f3Y-  -  -  pT,  (2.7) 

respectively,  for  the  spatial  and  temporal  dependence  of 
the  perturbations.  Here  the  separation  constant  u  is  an 
eigenvalue  to  be  obtained  through  solution  of  (2.6),  subject 
to  the  requirement  that  (2.3)  vanish  at  r  -  r^  and  r  = 

Once  v  is  known,  the  behavior  of  T(t)  is  completely 
determined  from  (2.6).  In  the  next  section  we  show  how  to 
solve  this  equation. 
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3.  Time  dependence  of  the  perturbations 

At  early  times,  when  f(t)  »1,  the  solutions  of  (2.7) 
are  approximately  given  by  exp[+  ivi*’t/T)J.  For  p<0,  one 
root  is  growing  and  one  is  damping,  while  for  ii>0,  both 
are  oscillatory.  Consequent!^  we  cap  expect  that  in  some 
sense  positive  u  will  be  associated  with  stable  behavior  and 
negative  p  with  instability. 

It  is  convenient  to  make  use  of  the  concepts  ot 
relative  amplification  and  relative  stability,  introduced 
by  BERNSTEIN  and  BOOK  (1978).  Suppose  that,  as  the 
imploding  pellet  or  shell  contracts,  the  size  of  the 
pertubations  decreases  also,  but  at  a  slower  rate.  That  is, 
f(t)  and  T(t)  both  approach  zero,  but  the  ratio  T(t)/f(t) 
grows,  thus  making  the  perturbations  more  desymmetrizing 
or  disruptive.  When  this  happens,  we  say  the  mode  is 
relatively  unstable.  Its  amplitude  can  eventually  exceed 
the  thickness  of  the  pellet  (by  which  point  the  linear 
approximation  is  no  longer  valid) . 

The  idea  of  "relative  stability"  should  be  applied  with 
a  degree  of  caution.  If  an  imploding  shell  is  ’perturbed" 
by  shifting  its  position  infinitesimally,  say  in  the 
direction  of  the  z  axis,  the  contraction  takes  place  just 
as  before.  A  linear  treatment,  however,  predicts  that  the 
i  m  3  0  mode  is  overstable  (i.e.,  it  grows  in  time 

without  oscillation) ,  evidently  a  spurious  conclusion.  On 
the  other  hand,  a  growing  oscillatory  mode  poses  a  genuine 
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threat  to  the  integrity  of  the  implosion,  even  if  the 
growth  is  only  relative  to  the  unperturbed  dimensions. 


When  VI  «•  1 ,  T  =  f  is  a  solution  of  (2.7)  .  The  other 

solution  is  found  by  the  method  of  variation  of  constants  to 

t 


be 

g (t)  =  T 

We  will  see  shortly  that  v>  =  l  is  the  critical  value  for 
determining  marginal  relative  stability. 

When  Y=l»  we  rewrite  (2.6)  in  terms  of  the  new 
independent  variable  x  =  i/i± . 


(3.1) 


d^T 

X  - r  + 


dx 


,1  .  dT  VI  „  - 

2  •  d7  -  2 


(3.2) 


a  confluent  hypergeometric  equation.  The  solutions  can 
be  expressed  as  a  linear  combination  of 

Tj  =  <KjVi;  j;  (3.3a) 

and 

T^  =  -(-2£/if)‘’  «(j  +  ju;  |-;£«f),  (3.3b) 

where  4(a;b;x)  is  the  Kummer  function  (see  e.g.,  hBRAMOWITZ 
and  STEGUN,  1964).  As  functions  of  t,  the  solutions  (3.3) 
satisfy  the  initial  condition 

Tj^(O)  -  1,  tTj^(O)  =  0>  (3.4a) 

T2(0)  -  0,  tT2(0)  =  1.  (3.4b) 
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for  large  negative  values  of  the  argument  M  f  both  behave 
as 


If  Y  “  1/  we  write  x  ■  1  -  f”“'  where  a  -  3(y-1).  Then 
(2.6)  becomes  the  hypergeometric  equation 


<j^T 

x(l-x)  + 

dx 


['  -  (•♦‘•♦Dx  I;  - 


abT  «  0 . 


(3.5) 


Here 


a 

b 


(r±A)/4o, 


(3.6a,b) 


2  ^ 

where  p  =  “  +  2  and  A  ■  (T  -  Bay)  ,  and  c  =  >s .  The 


solutions  of  (3.5)  can  be  chosen  as 


T^  -  F(a,b; H  ;  l-f"“)  ,  (3.7a) 

T2  =  -  f»  l-f’“)  ,  (3.7b) 

where  F(a,b;c;x)  is  Gauss'  hypergeometric  function.  The 

solutions  (3.7)  again  satisfy  (3.4).  As  f  0,  both  approach 

asymptotic  forms  (ABRAMOWITZ  and  STEGUN,  1964)  containing 

terms  proportional  to  f  with  the  lower  sign,  this 

expression  diverges  as  f  -»•  0  whenever  A  >  T,  i.e., 

y  <  0.  Since  the  condition  for  A  to  be  real  is  that  y 

2 

be  no  greater  than  6(a+2)  /8a;  whose  minimum  value 
as  a  function  of  a  is  unity,  we  see  that  y  <  1  is  always 
sufficient  to  malce  T/f  diverge.  If  A  is  imaginary,  T/f 
still  diverges  when  F  <  4  (y  <  5/3).  BOOK  (1978)  has 
presented  a  simple  argument  involving  conservation  of  wave 
action  to  explain  this  result  in  terms  of  the  geometrical 
properties  of  the  implosion. 
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4.  Rayleigh-Taylor  instability 


He  turn  now  to  the  problem  of  calculating  p.  He  dot 
(2.5)  with  ^  and  use  (1.6)  to  rewrite  the  coefficients 
(except  in  the  first  term)  in  terms  of  p  Instead  of  p. 
Introducing  the  notation  a=V and  X  and  using  the 


identities 


•  Vp  =  7  •  (po^)  -  po'^  -  p^  •  7o,  (4.1) 

5  •  7(po)  =  7  •  (po^)  -  po^,  (4.2) 

^  •  (70  •  7p  =  7  •  (pi^  •  7^)  +  -pu»^  (4.3) 

where  the  transpose  is  defined  in  Cartesian  coordinates  by 
ij  = 

T  ^pp€^  =  P  [(y  -  1)  +  7^  s  -  u>^]+  yV  •  (po^).  (4.4) 

Integration  over  the  volume  occupied  by  the  shell  yields 


pj*  dVp{;^  =  j  dVp  j^(Y-l)o^  +  7^ 
V  V 

+  Y  /  dS  pc^  .  ^ 


: 


(4.5) 


where  in  the  last  term,  the  unit  vector  n  is  so  defined  that 
it  points  away  from  the  shell  on  both  inner  and  outer 
surfaces.  The  left  hand  member  and  the  first  terms  in  the 
volume  integral  on  the  right  hand  side  of  (4.5)  are  manifestly 
nonnegative.  From  this  we  see  that  relative  instability 
(p  <  1)  can  result  in  only  two  ways:  if  7  x^  ¥  0,  or  if 
the  integrated  terms  are  nonvanishing.  The  latter  ii?  the 


a;..  !,,.....  _ 


case  whenever  a  perturbation  exists  at  a  point  where  the 
density  changes  discontinuously .  For  the  model  we  have 
assumed  in  S  2,  one  such  point  is  always  located  at  the 
shell  outer  radius,  r  «  "^0^*** 

The  external  pressure,  which  enters  the  model  as  a 
boundary  condition,  produces  an  inward  acceleration.  There 
is  thus  an  effective  gravity  in  the  outward  direction. 
Hence  we  anticipate  that  a  Rayleigh-Taylor  instability 
should  occur,  localized  at  the  outer  surface. 

In  slab  geometry  under  uniform  destabilizing 
gravitational  acceleration  G,  the  Rayleigh-Taylor  growth 
time  for  incompressible  modes  is  given  by 


-2 


vr  "*2-^1  , 

kG  — T -  ' 


(4.6) 


where  k  is  the  wavenumber  and  p  and  p  are  the  densities 

2  1 

above  and  below  the  Interface.  In  the  present  problem, 

pj^  »  0  and  we  can  write  k  =  2ir/X,  where  the  wavelength  satisfies 

X  =  2irr/f.  A  reasonable  approximation  to  the  effective 
-2 

gravity  is  G  ^rx  .  Substituting  in  (4.6),  we  thus  estimate 


that  the  growth  time  will  be  given  by 


-  T^/i.  (4.7) 

9 

To  study  the  analogous  modes  in  our  model  of  a 
uniformly  imploding  shell,  we  begin  by  operating  on  (2.5) 
with  the  divergence  and  with  the  curl.  The  results  are 
given  by 
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■  4l--^ 


7 


ua  +  y{t^[p  ^  7^(po)  +  V  (p”^)  *7  (pa)J 


-2a  -  ^  •  V  M  ig}  0 


(4.8) 


(y-1)  u  +^X7o  +  t^y7(p  X7(pa)  -  0.  (4.9) 

From  (4.8)  -  (4.9)  it  follows  that  ^  vanishes  if  o  does. 


If  that  happens. 


^ 


(4.10) 


where  the  potential  ((i  satisfies  Laplace's  equation. 


7  (fr  «  0 . 


(4.11) 


The  general  solution  of  (4.11)  In  spherical  coordinates  Is 
a  sum  of  terms  of  the  form 

<(i(r,0,c(l  *  ((|i^r*'  +  <l>_r"^"^)  Y^^(0  ,cp)  ,  (4.12) 

where  are  constants.  Now  from  (4.10)  and  (4.11),  (2.6) 

becomes 

(y-l)j^  +7(jc*y  =  7  j^(y-l)  (J>  +  ^  .  Vi|.J  *  0.  (4.13) 

Hence  either  y  =  *-+2  or  p  =  -  f+l,  corresponding  to  the 
two  terms  In  (4.12).  The  latter  mode,  localized  (not 
surprisingly)  at  the  outer  surface  of  the  shell.  Is 
unstable.  From  (2.7)  we  see  that  the  time  scale  for  this 


mode  is 


Tg  -  T^/d+l)  , 


(4.14) 


i 


essentially  the  same  as  the  estimate  given  in  (4.7)  . 

Note  that  this  result*  Identical  with  that  obtained 
by  KIDDER  (1976)  for  the  case  of  isentropic  implosions 
with  Y  *  5/3, is  completely  independent  of  any  details  of 
the  density,  pressure  or  entropy  profile,  the  magnitude 
of  r^  or  d,  and  the  value  of  y.  The  unstable  mode  is 
purely  Incompressible,  even  though  the  unperturbed 
fluid  is  experiencing  compression  under  the  influence  of 
the  implosive  force. 
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5.  Convective  Instability 

Ne  consider  perturbations  for  which  a  does  not 
vanish.  In  that  case  we  see  from  (4.9)  that  u  is  also 
nonzero,  unless  a  depends  only  on  r  or 


r  +  T  YP^(P  )  =  0,  (5.1) 


which  is  equivalent  to 


/  “Y. 

57  -  <pp  )  -  0. 


(5.2) 


The  former  possibility,  which  amounts  to  looking  just  at 
i-0  modes,  and  the  latter,  which  was  discussed  by  BOOK 
(1978) ,  are  both  included  in  the  general  treatment  to  be 
presented  here. 

We  begin  by  moting  that  ja  enters  (4.8)  only  in  the 
form  v=  An  expression  for  this  quantity  can  be  found 

by  taking  the  curl  of  (4.9),  then  dotting  with^.  Eliminating 

V  in  favor  of  o  and  writing  o(r)  =  a  (r)  t,  (0,0)),  we 

ilm 

obtain  an  ordinary  differential  equation 

T  ^(vi-2)o  +  YP  ^  ^po)  +  2r"^(po)  -  i  ( i  +  1 )  r“  ^poj 

+  Y(P"^)’  (po)'  -  '  <5.3) 

where  the  prime  denotes  differentiation  with  respect  to  r. 
Equation  (5.3)  can  be  rewritten  in  the  form 

where  w  >  -  po  is  the  radial  part  of  p ^  and 


0  (5.4) 
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^  -  f 


(5.5) 


After  multiplying  (5.5)  by  w  and  integrating  with 
respect  to  r,  we  find 


I  dr  r^p  ^w'^  ~  J  dr  r^p  ^[1(1+1) 
rn 


i-k-  P_L  .  i  )  +  KEl]  w"  .  0 

%-l  ro  2’  2rp-' 


.  _  2 


(5.6) 


‘y-1  rp  ^2'  2rp’ 


Here  the  surface  terms  left  after  the  integration  by  parts 
vanish  because  the  perturbed  pressure  must  be  zero  at  the 
boundaries.  This  equation  is  self-adjoint  and  almost  in 
Stura*Liouville  form.  Extremization  subject  to  the 
requirement  that  p  be  stationary  recovers  the  differential 
equation  (5.4),  so  that  (5.6)  is  the  basis  for  a  varia¬ 
tional  principle. 


Defining 


dr  r  w  p*/pp» 


(5.7a) 


t  .  -  »  *  yjir  +  t(»,+  l)/ar  w^/p»  (5.7b) 

Cm.  A(i  +  1)  (A  +  yJax  X  w^pVP^)  (5.7c) 


(5.7c) 


where  all  integrals  run  from  to  tQ+d,  we  can  express 
(5.6)  as 


The  eigenvalue  y  accordingly  satisfies 


y-1 


-B  ±  (B^-4AC)** 
2A 


(5.9) 


From  the  definitions,  B  >  0  >  A.  Thus  the  lower  root  always 
satisfies  y  >  1,  ruling  out  all  but  the  weak  "geometric" 
type  of  amplification  (§3).  We  can  make  B  arbitrarily 
large  without  changing  A  or  C  significantly  by  superposing 
a  rapid  Infinitesimal  "wiggle"  on  w.  As  B  Increases 
without  bound,  the  upper  root  approaches  y«  1  .  Hence  we 
see  that  the  extremal  y  corresponding  to  the  exact 
eigenfunction  w  Is  a  minimum .  It  follows  that  there  Is  an 
unstable  mode  (y  <  1)  If  any  test  function  w  can  be  found 
for  which  (5.9)  yields  y  <  1.  The  condition  for  this  to 
occur  Is  C  >  0.  Since  the  latter  can  be  written 


C  ■  l(i+l) y  dr  r  w^  s'/sp,  (5.10) 

a  necessary  and  sufficient  condition  that  y  <  1  Is  s'(r)>  0 

somewhere . 

The  mechanism  for  this  Instability  was  described  by 
BOOK  (1979) .  Because  the  unperturbed  system  Is  nonsteady 
(f  f  0),  individual  elements  of  fluid  are  subject  to  an 
effective  gravitational  acceleration  G  ■  -rf,  so  that 
elements  of  differing  density  experience  different  buoyancy 
forces.  If  the  result  of  Interchanging  two  differential 


elements  with  different  initial  radii  is  to  reduce  the  sum 
of  the  kinetic,  compresslonal,  and  effective  gravitational 
energies  associated  with  these  elements,  the  system  is 
unstable.  It  can  easily  be  shown,  by  modifying  slightly 
the  argument  of  BOOK  (1979),  that  the  condition  for  this  to 
occur  is  just  that  the  profiles  be  " super adiabatic , "  i.e., 
that  s'(r)  >  0  somewhere  in  the  shell.  We  thus  anticipate 
that  the  instability  associated  with  this  condition  will  be 
characterized  by  overturning  of  profiles,  such  as  is 
typically  seen  in  convective  instabilities  of  static  media 
(LANDAU  and  LIFSHITZ,  1959).  There  also  the  condition  for 
instability  is  the  occurrence  of  superadiabaticlty . 

Note  that  internal  Rayleigh-Taylor  modes,  i.e,  modes 
associated  with  the  occurrence  of  opposing  pressure  and 
density  gradients  (p'p  '  <  0)  at  some  point  r  satisfying 
r^  <  r  <  r^  ■t'd,  are  described  by  this  formalism,  since  now 
both  terms  in 

s'/a  -  p'/p  -  YP’/P  (5.11) 


are  positive 
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where  p  is  constant  and  L  is  a  characteristic  scale  size. 
If  p(0)  ■  0  it  follows  that 


-  .  .  -  2  , ,  2  ,  2, 

p(r)  =  pr  /{r  +L  ) 


(6.11) 


s(r)  *  p  p  (r  /L  ) (1+r  /L  ) 


(6.12) 


so  that  again  s’(r)  >0.  The  logarithmic  derivatives  of  p 

and  p  are 


P.  '  -4r/L'^ 

^  “  l.r^/L^ 


(6.13) 


^  r(l  +  r^/L^) 

If  we  let  w(r)  =  r^y(r),  where 


(6.14) 


q  =  +  [<2i  +  l)^  -  4Kj^|,  (6.15) 


then  y  satisfies  the  hypergeometric  equation  (3.5)  with 


-4q;j  (6.16) 


2  2 

c  »  q  +  3/2,  and  x  =  -  r  /L  .  The  condition  that  w(r)  vanish 
at  infinity  reduces  to  q  -  2b  <  0 ,  or  p  >  -3,  independent  of 
i.  The  discreet  set  of  unstable  modes  has  thus  merged  into 
a  continum,  the  result  of  moving  the  boundary  to  infinity. 
Finally,  letting  =  0  but  keeping  d  finite,  choose 


/  2,,  2. 

p  =  p  exp  (r  /L  ) 


(6.17) 


T  hen 


P  =  p  exp  (r^/L^) , 


(6.18) 
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and 


£'=  iL  .  (6.19) 

P  P  ^2 

For  this  choice,  s(r)  is  monotone  decreasing ; 

s(r)  =  pp~Y  exp  ^(1-Y)r^/L^|  .  (6.20) 

I 

Writing  w(r)  ^r  z(r),  we  find  that  z  satisfies  the 
confluent  hypergeometric  equation 


xz"  +  (c-x)  z'  -  az  =  0  , 


(6.21) 


where 


a 


1  iLzl  , 

2  V-1 


(6.22) 


2  2 

c  «  £'f3/2,  and  x  *  r  /L  .  From  the  series  expansion,  it  is 

2  2 

clear  that  the  Kummer  function  t(a;c;x)  can  vanish  at  x  «  d  /L 
only  If  a  <  0.  It  follows  that 

W-1  >  h  {yl*l  ±J(y£  +  1)^  -  41(1  +  1)  (Y-1)]  **  )  >  0,  (6.23) 

so  that  all  modes  are  stable,  as  may  be  inferred  from  the 
fact  that  s' (r)  <  0. 
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Discussion 


( 


i: 


7  . 

We  have  seen  that  for  a  broad  class  of  pressure  and 
density  profiles  it  is  possible  to  calculate  the 
evolution  of  linear  perturbations  on  a  uniformly  imploding 
pellet  or  shell.  The  calculations  presented  above  show 
that  the  time  development  of  the  perturbations  depends  only 
on  the  eigenvalue  ii  and  y,  the  ratio  of  specific  heats.  If 
p  >  1,  the  perturbation  amplitude  grows  weakly  relative  to 
the  radius  of  the  unperturbed  pellet  during  implosion  if  and 
only  if  Y  <  5/3.  If  p  <  1,  perturbation  amplitudes  are 
nonoscillatory  and  are  always  relatively  amplified,  the 
degree  of  amplification  decreasing  with  y.  When  p  <  o  the 
mode  amplitudes  grow  absolutely , not  just  relatively.  They 
grow  rapidly  if  the  eigenvalue  p  is  large  and  negative.  These 
results  are  summarized  in  Table  1,  which  presents  the  result 
of  integrating  (2.7)  for  various  choices  of  y  and  p  until 
f  =  0.1  (tenfold  radial  compression) , and  in  Table  2,  which 
shows  the  corresponding  results  for  f  =  0.05  (twentyfold 
compression) . 

Incompressible  irrotational  modes,  which  are 
physically  of  Ray leigh-Taylor  type,  have  p  =  1  ~  t , 
irrespective  of  the  shape  of  density  and  pressure  profiles. 

For  compressible  (convective)  modes,  p  <  1  or  p  <  0  can  also 
occur,  but  the  limiting  value  as  t  is  finite.  These 

modes  are  consequently  less  significant  unless  (i)  the 
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compressible  disturbances  are  Initiated  with  larger 
amplitude}  or  (11)  some  change  In  the  physical  situation 
introduces  a  stabilizing  influence  at  the  outer  boundary, 
where  the  Ray leigh-Tay lor  modes  are  localized,  but  does  not 
affect  the  convectively  unstable  region  where  s’(r)  >  0. 

For  real  shells,  which  are  affected  by  ablation. 

Ionization,  thermal  conduction  and  other  processes  absent 
from  our  model,  these  results  can  only  be  applied  in  a 
rough  approximation.  The  model  Is  most  nearly  valid  in 
the  interior  of  the  shell,  far  from  the  deposition  zone. 

In  this  region  it  should  predict  correctly  the  stability  of 
the  target  as  implosion  and  compression  take  place. 

i 

A  reasonable  approximation  is  to  take  y  4/3  and 
Inquire  regarding  the  effect  of  a  radial  compression  of 
10  to  20,  corresponding  to  a  density  Increase  of  three 
to  four  orders  of  magnitude.  As  p  decreases,  the 
oscillation  period  and  growth  rate  both  Increase.  From 
numerical  Integration  of  (2.6)  for  various  choices  of  p 
between  -3  and  3  (Tables  1  and  2)  It  Is  clear  that  relative 
amplification  by  a  factor  of  100  Is  quite  possible.  This 
Imposes  a  limitation  on  acceptable  asymmetries  or 
Irregular ities , restricting  them  to  relative  magnitudes  of 
order  1%  or  less.  Even  the  convective  modes,  for  which 
-p  Is  bounded  above,  can  be  amplified  significantly  relative 
to  the  Imploding  shell  raduls.  Of  course  the  Raylelgh- 
Taylor  modes  at  large  values  of  t  are  amplified  much  more. 
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since  only  the  simplest  properties  of  the  target 
material  were  Invoked  in  analysing  stability  of  the 
Implosion,  the  predicted  behavior  should  be  observed 
universally.  That  is,  it  does  not  depend  on  details 
of  the  deposition  process  or  on  subtle  plasma  mechanisms. 
(Of  course,  if  the  target  becomes  ionized  in  the  course  of 
implosion,  the  model  must  be  altered  to  reflect  the 
ensuing  complications.)  Viewed  in  this  light,  the 
present  calculation  provides  a  zeroth  order  prediction  of 
the  destabilizing  and  desynmetr izing  effects  of  implosion 
of  a  compressible  target. 
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